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Abstract 

We study diffuse gamma-ray emission at intermediate Galactic latitudes measured by the Fermi Large Area Tele- 
scope with the aim of searching for a signal from dark matter annihilation or decay. In the absence of a robust 
dark matter signal, we set conservative dark matter limits requiring that the dark matter signal does not exceed the 
observed diffuse gamma-ray emission. A second set of more stringent limits is derived based on modeling the fore- 
ground astrophysical diffuse emission. Uncertainties in the height of the diffusive cosmic-ray halo, the distribution of 
the cosmic -ray sources in the Galaxy, the cosmic -ray electron index of the injection spectrum and the column density 
of the interstellar gas are taken into account using a profile likelihood formalism, while the parameters governing the 
cosmic -ray propagation have been derived from fits to local cosmic -ray data. The resulting limits impact the range 
of particle masses over which dark matter thermal production in the early Universe is possible, and challenge the 
interpretation of the PAMELA/Fermi-LAT cosmic ray anomalies as annihilation of dark matter. 

Keywords: 



1. Introduction 

The Milky Way halo has long been considered a good 
target for searches of indirect signatures of dark matter 
(DM). WIMP DM candidates are expected to produce 
gamma rays, electrons and protons in their annihilation 
and decays and such emission originating in our Galaxy 
would appear as a diffuse signal. At the same time, the 
majority of the Galactic diffuse emission is produced 
through radiative losses of cosmic -ray (CR) electrons 
and nucleons in the interstellar medium. Modeling of 
this emission presents one of the major challenges when 
looking for subdominant signals from dark matter. 



In this analysis we test the diffuse LAT data for a con- 
tribution from the DM signal by performing a fit of the 
spectral and spatial distributions of the expected pho- 
tons at intermediate Galactic latitudes. In doing so, we 
take into account the most up-to-date modeling of the 
established astrophysical signal, adapting it to the prob- 
lem in question (TJ. Our aim is to constrain the DM 
properties and treat the parameters of the astrophysical 
diffuse gamma-ray background as nuisance parameters. 
Besides this approach, we will also quote conservative 
upper limits using the data only (i.e. without performing 
any modeling of the astrophysical background). 
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We follow [2| in using the GALPROP code v54 0, 
to calculate the propagation and distribution of CRs in 
the Galaxy and the whole sky diffuse emission. In 
various standard parameters of the CR propagation were 
studied in a fit to CR data and it was shown that they rep- 
resent well the gamma-ray sky, although various resid- 
uals (at a ~ 30% level [2|), both at small and large 
scales, remain. In our work, we use the results of the 
fits to the CR data from [ 2 1 but we allow for more free- 
dom in certain parameters governing the CR distribu- 
tion and known astrophysical diffuse emission and con- 
strain these parameters by fitting the models to the LAT 
gamma-ray data. Despite the large freedom we leave 
in the models we see residuals in our ROI at the +30% 
level and at ~ 3 cr significance. These residuals can 
be ascribed to various limitations of the models: imper- 
fections in the modeling of gas and ISRF components, 
simplified assumptions in the propagation set-up, un- 
resolved point sources, and large scale structures like 
Loop I [4 1 or the Galactic Bubbles [5 |. Since residuals 
do not seem obviously related to DM, we focus in the 
following on setting limits on the possible DM signal, 
rather than searching for a DM signal. 

2. DM maps 

We parametrize the smooth DM density p with a 
NFW spatial profile Q 



PoR s 



r (1 + r/R s ) 2 
and a cored (isothermal-sphere) profile 



P(r) 



{r 2 +RD ' 



(1) 



(2) 



For the local density of DM we take the value of po = 
0.43 GeV cm" 3 (SJ, and the scale radius of R s = 20 kpc 
(for NFW) and R c = 2.8 kpc (isothermal profile). We 
also set the distance of the solar system from the cen- 
ter of the Galaxy to the value R G = 8.5 kpc. For the 
annihilation/decay spectra we consider three channels 
with distinctly different signatures: annihilation/decay 
into the bb channel, into p + p~, and into r + r~. In the first 
case gamma rays are produced through hadronization of 
annihilation products and subsequent pion decay. The 
resulting spectra are similar for all channels in which 
DM produces heavy quarks and gauge bosons and this 
channel is therefore representative for a large set of par- 
ticle physics models. The choice of leptonic channels 
provided by the second and third scenarios, is motivated 
by the dark matter interpretation [9| of the PAMELA 



positron fraction ifTUll and the Fermi LAT electrons plus 
positrons ifTTI measurements. In this case, gamma rays 
are dominantly produced through radiative processes of 
electrons, as well as through the Final State Radiation 
(FSR). We produce the DM maps with a version of 
GALPROP slightly modified to implement custom DM 
profiles and injection spectra (which are calculated by 
using the PPPC4DMID tool described in ifTH '). 

3. Approach to set DM limits 

We use 24 months of LAT data in the energy range 
between 1 and 100 GeV (but, we use energies up to 
400 GeV when deriving DM limits with no assumption 
on the astrophysical background). We use only events 
classified as gamma rays in the P7CLEAN event selec- 
tion and the corresponding P7CLEAN_V6 instrument re- 
sponse functions (IRFsJ] Structures like Loop I and 
the Galactic Bubbles appear mainly at high Galactic lat- 
itudes and to limit their effects on the fitting we will 
consider a ROI in Galactic latitude, b, of 5° < \b\ < 15°, 
and Galactic longitude, /, |/| < 80°. We mask the region 
\b\ < 5° along the Galactic Plane, in order to reduce the 
uncertainty due to the modeling of the astrophysical and 
DM emission profiles. 

3.1. DM limits with no assumption on the astrophysical 
background 

To set these type of limits we first convolve a given 
DM model with the Fermi IRFs to obtain the counts ex- 
pected from DM annihilation. The expected counts are 
then compared with the observed counts in our ROI and 
the upper limit is set to the minimum DM normalization 
which gives counts in excess of the observed ones in at 
least one bin, i.e. we set 3cr upper limits given by the 
requirement n iDM - 3 yjn iDM > n,-, where n iDM is the ex- 
pected number of counts from DM in the bin i and «, the 
actual observed number of counts. 

3.2. DM limits with modeling of astrophysical back- 
ground 

In this analysis we model the diffuse emission as a 
combination of a dark matter and a parameterized con- 
ventional astrophysical signal and we derive the lim- 
its on the DM contribution using the profile likelihood 
method. More precisely, for each DM channel and mass 
the model which describes the LAT data best is the one 



http://fermi.gsfc.nasa.gov/ssc/ 
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which maximizes the likelihood function defined as a 
product running over all spatial and spectral bins i , 

Lic(0dm) = L k (0DM, a) = maxg ]~ [ P ik {nr, a, 8 DM ), (3) 

i 

where P,k is the Poisson distribution for observing «, 
events in bin i given an expectation value that depends 
on the parameter set (9dm, <?)■ 8dm is the intensity of 
the DM component, a represents the set of parameters 
which enter the astrophysical diffuse emission model as 
linear pre-factors to the individual model components 
(cf. equation [4] below), while k denotes the set of pa- 
rameters which enter in a non-linear way. We sample 
non-linear parameters of the astrophysical background 
on a grid (for computational efficiency). 

The linear part of the fit is performed with the 
GaRDiAn code in the following way. The CR source 
distribution (CRSDs) is a critical parameter for DM 
searches and we define a parametric CRSD as sum of 
step functions in Galactocentric radius R, and treat the 
normalization of each step as a free parameter in the fit 
to gamma ray^ In order to have conservative and ro- 
bust limits we set to zero the e, p CRSDs in the inner 
Galaxy region, within 3 kpc of the Galactic Center. In 
this way, potential e and p CR sources which would be 
required in the inner Galaxy will be potentially compen- 
sated by DM, producing conservative constraints. The 
linear part of the fit is then performed by combining 
the sky maps produced in each ring of CRSD of the 
main components of the diffuse emission: the emission 
from 7T° decay, bremsstrahlung and inverse Compton. 
Additionally an isotropic component arising from the 
extragalactic gamma-ray background and misclassified 
charged particles needs to be included to fit the Fermi 
LAT data. In summary, the various GALPROP outputs 
are combined as: 

j ' 
a x (Xy + Xic) + aiGB - m IGB "' ■ (4) 

m 

The sum over i is the sum over all step-like CRSD 
functions, the sum over j corresponds to the sum over 
all Galactocentric annuli (details of the procedure of a 
placement of the gas in Galactocentric annuli and their 



-CRSDs are traditionally modeled from the direct observation of 
tracers of SNR and can be observationally biased. 



boundaries are given in J2JI) • H denotes the gamma-ray 
emission from atomic and ionized interstellar gas while 
H2 the one from molecular hydrogerj^] and IC the In- 
verse Compton emission. \y w &Xic are the prompt and 
Inverse Compton (when present) DM contribution and 
a x the overall DM normalization. aicB,m and IGB m de- 
note the Isotropic Gamma-ray Background (IGB) inten- 
sity for each of the five energy bins over which the index 
m runs. In all the rest of the expression the energy index 
m is implicit since we don't allow for the freedom of 
varying the GALPROP output from energy bin to energy 
bin. We do not include sources in the fit as we use a 
mask to filter point sources from the 1FGL catalog [13]. 

The outlined procedure is then repeated for each set 
of values of the non-linear propagation and injection 
parameters to obtain the full set of profile likelihood 
curves. We scan over the three parameters: electron in- 
jection index, the height of the diffusive halo and the 
gas to dust ratio which parametrizes different gas col- 
umn densities (see Table[T]i. In this way we end up with 
a set of k profiles of likelihood L^Adm), one for each 
combination of the non-linear parameters. The envelope 
of these curves then approximates the final profile likeli- 
hood curve, L(6dm), where all the parameters, linear and 
non-linear have been included in the profile. Limits are 
calculated from the profile likelihood function by find- 
ing the DM ,iim values for which L(6 DM ,iim)/L(0DM,max) is 
exp(-9/2) and exp(-25/2), for 3 and 5 <x C.L. limits, 
respectively. 

4. Results 

An important point to note is that, for each DM 
model, the global minimum we found lies within the 
3(5) <x regions of many different models, providing a 
check against a bias in our procedure. This point is il- 
lustrated in Figure [Tj where the profile likelihoods for 
the three nonlinear parameters, Zh, Je,z an d d2HI, are 
shown. To ease reading of the figure the profiling is 
actually performed with further grouping DM models 
with different DM masses, but keeping the different 
DM channels, DM profiles and the annihilation/decay 
cases separately. The curve for the fit without DM is 
also shown for comparison. Each resulting curve has 
been further rescaled to a common minimum, since we 
are interested in showing that several models are within 
-2AlogL < 25 around the minimum for each DM fit. 



3 It should be noted that in our case, where we mask ±5° along the 
plane, the expression actually simplifies considerably since only the 
local ring Xqo factor enters the sum, since all the other Hi rings do 
not extend further than 5 degrees from the plane. 
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The y e ,2 profile, for example, indicates that all mod- 
els with y Cj 2 from 1.9 to 2.4 are within -2AlogL < 
25 around the minimum illustrating that the sampling 
around each of the minima for the six DM models is 
dense. Similarly, the d2HI profile indicates that all mod- 
els with d2HI in the range (0.120 - 0.160) xl0~ 20 mag 
cm 2 are within 5<x from the minima for each of the six 
DM models. Finally the Zh profile indicates that basi- 
cally all the considered values of Zh are close to the ab- 
solute minima. This last result is not surprising since, 
within our low-latitude ROI, we have little sensitivity to 
different Zh and basically all of them fit equally well. 
There is some tendency to favor higher values of Zh 
when DM is not included in the fit, while with DM the 
trend is inverted. Although the feature is not extremely 
significant it is potentially very interesting. 

Upper limits on the velocity averaged annihilation 
cross section into various channels are shown in Fig. [2] 
for isothermal profile of the DM halo^] together with 
regions of parameter space which provide a good fit 
to PAMELA (purple) and Fermi LAT (blue) CR elec- 
tron and positron data [14|. Despite the various conser- 
vative choices described above, the resulting DM lim- 
its are quite stringent. They are comparable with the 
limits from LAT searches for a signal from DM an- 
nihilation/decay in dwarf galaxies 1151 . In particular, 
as shown in Figure [2] for masses around 20 GeV the 
thermal relic value of the annihilation cross section is 
reached, both for the bb and t + t~ channels. 

Overall, rather than being due to residual astrophysi- 
cal model uncertainties, the remaining major uncertain- 
ties in the DM constraints from the Halo region come 
from the modeling of the DM signal itself. The main 
uncertainty is in the normalization of the DM profile, 
which is fixed through the local value of the DM den- 
sity. We use the recent determination p = 0.43 GeV 
cirT 3 from |8|, which has, however, a large uncertainty, 
with values in the range 0.2-0.7 GeV cm -3 still viable. 
A large uncertainty in po is particular important for an- 
nihilation constraints since they scale like pi, while for 
constraints on decaying DM the scaling is only linear. 
A less important role is played by the uncertainties in 
the DM profile, since in our region of interest different 
profiles predict similar DM densities. A better determi- 
nation of the local DM density, as well as of the param- 
eters determining the global structure of the DM Halo, 
is therefore of the utmost importance for reducing the 
uncertainties related to DM constraints from DM halo, 
but it is beyond the scope of this paper and is the subject 
of dedicated studies. 
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Limits obtained using the NFW profile are only slightly better. 
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Non linear Parameters 


Symbol 


Grid values 


index of the injection CRE spectrum 
half height of the diffusive halcj^j 
dust to HI ratio 


Zh 

d2HI 


1.800, 1.925, 2.050, 2.175, 2.300, 2.425, 2.550, 2.675 
2, 4, 6, 8, 10, 15 kpc 
(0.0120, 0.0130, 0.0140, 0.0150, 0.0160, 0.0170) xl0~ 20 mag cm 2 


Linear Parameters 


Symbol 


Range of variation 


eCRSD and pCRSD coefficients 
local H 2 to CO factor 
IGB normalization in various energy bins 
DM normalization 


e P 
C C 

vloc 
A CO 

OlCB,m 

a x 


0,+oo 

0-30 xlO 20 cm" 2 (K km s" 1 )- 1 
free 
free 



Table 1 : Summary table of the parameters varied in the fit. The top part of the table show the non linear parameters and the grid values at which the 
likelihood is computed. The bottom part show the linear parameters and the range of variation allowed in the fit. The coefficients of the CRSDs are 
forced to be positive, apart and c^ which are set to zero. The local Xco ratio is restricted to vary in the range 0-30 XlO cm 
while aiGB.m and a x are left free to assume both positive and negative values. See the text for more details. 





Figure 1: Profile likelihood curves for zh, 7e,2 and d2HI. The various curves refer to the case of no DM or different DM models (see the legend in 
the figure, where we mark a dominant decay (DEC) or annihilation (AN) channel and the assumed DM profile). All minima are normalized to the 
same level. Horizontal dotted lines indicate a difference in -2AlogL from the minimum of 9 (3cr) and 25 (5<x). 
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Figure 2: Upper limits on the velocity averaged DM annihilation cross-section (left) and decay time (right) including a model of the astrophysical 
background compared with the limits obtained with no modeling of the background. Limits are shown for bb (upper), p^yT (middle) and t + t~ 
(lower panel) channels, for a DM distribution given by the isothermal distribution. The horizontal line marks the thermal decoupling cross section 
expected for a generic WIMP candidate. The regions of parameter space which provide a good fit to PAMELA (purple) and Fermi LAT (blue) CR 
electron and positron data are also shown. 



